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Abstract. We non-perturbatively analyze the effect of electron-electron interactions 
on weak localization (WL) in relatively short metallic conductors with a tunnel barrier. 
We demonstrate that the main effect of interactions is electron dephasing which persists 
down to T = and yields suppression of WL correction to conductance below its non- 
interacting value. Our results may account for recent observations of low temperature 
saturation of the electron decoherence time in quantum dots. 
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Electrons propagating in a disordered conductor get scattered and interfere. This 
quantum interference is possible only as long as the electron wave functions remain 
coherent. In any realistic situation, however, interactions between electrons and with 
other degrees of freedom may limit phase coherence and, hence, reduce electrons ability 
to interfere. The interplay between scattering, quantum coherence and interactions 
yields a rich variety of non-trivial effects and significantly impacts electron transport in 
disordered conductors. 

The so-called weak localization (WL) correction to the conductance of a disordered 
system Gwl is most sensitive to electron coherence and is known to arise from 
interference of pairs of time- reversed electron paths [1]. In a system of two scatterers 
separated by a cavity (quantum dot) and in the absence of interactions this correction 
can be directly evaluated [2J. The effect of electron-electron interactions can be described 
in terms of fluctuating voltages. Provided the voltage drops only across the barriers and 
not inside the cavity electron-electron interactions yield energy dependent logarithmic 
renormalization of the dot channel transmissions |3l H] but do not cause any dephasing 
[5|, E]. The latter result can easily be understood if one observes that the voltage- 
dependent random phase acquired by the electron wave function \E' along any path turns 
out to be the same as that for its time-reversed counterpart. Hence, in the product \E'\E'* 
these random phases cancel each other exactly and quantum coherence of electrons 
remains preserved. 

It is important, however, that this cancellation occurs only in the case of two 
scatterers, whereas in a system of three or more scatterers the situation is entirely 
different. Consider, e.g., a system of two quantum dots depicted in Fig. 1 and again 
assume that fluctuating voltages are concentrated at the barriers. The phase factor 
accumulated along the path (see Fig. 1) which crosses the central barrier twice (at 
times ti and t > ti) and returns to the initial point (at a time tf) is e*''^^'-**-'^'^^''*^], 
where <f~^/e = V{t) is the fluctuating voltage across the central barrier. Similarly, the 
phase factor picked up along the time-reversed path reads e*''^^^*-''"'"**"*-'"'^^^*^'". Hence, 
the overall phase factor acquired by the product for a pair of time-reversed paths 
is exp(i$tot), where ^tot{ti,tf,t) = (p^iti) - ip+{t) - (p+{tf + ti-t) +(p+(tf). Averaging 
over phase fluctuations, which for simplicity are assumed Gaussian, we obtain 

^ ^-2F{t~ti)-2F{tf-t)+F(t;-U)+F(t;+U-2t) ^^-^ 

where we defined the phase correlation function 

F(t) = ((^+(t)-y,+ (0))^)/2. (2) 

Should this function grow with time the electron phase coherence decays and Gwl gets 
suppressed below its non-interacting value. 

The above arguments are not specific to the system of three barriers and after proper 
generalization can be applied to virtually any disordered conductor. At the same time, 
these arguments are not yet sufficient to quantitatively describe the decoherence effect 




of electron-electron interactions for two important reasons: (i) fluctuating voltages are 
treated as external (classical) fields rather than quantum fields produced internally by 
fluctuating electrons and (ii) Fermi statistics is not yet accounted for. Below we will 
cure both these problems and non-perturbatively evaluate WL correction Gwl for a 
metallic system with a tunnel barrier and (at least) two more scatterers in the presence 
of electron-electron interactions which turn out to reduce phase coherence of electrons 
at any temperature down to T = 0. 

We will consider a system with a tunnel barrier with dimensionless conductance 
Qt, which separates two sufficiently short disordered metallic conductors with Thouless 
energies E'xh and dimensionless conductances gL,R ^ ^^Qt- This system is described by 
the Hamiltonian 

H = HL + HR + f + H,^, (3) 

where Hl,r = 'Ea=hl lL,Rd^f '^L,L,Ri'^)HL,R^a;L,R{'^) IS the Hamihonian of the left 
(right) lead, Hl,r = ~ 2mL r ~ ~^ Uimp{r) is the single electron Hamiltonian in the 



left (right) lead, and T = X)a=t i Ij ^^'^ 



t(r)e"*'^(*)^^(r)^„;ij(r) + c.c. is the tunnel 
Hamiltonian. Here is the phase operator, which is related to the voltage drop across 
the junction 0(t) = J^^dt'eV{t'), and the r-integration runs over the junction area. 
Finally, Hem oc V'^ is the quadratic Hamiltonian of electromagnetic fields, the precise 
form of which depends on the circuit configuration and will not be specified here. 

Following the standard procedure we integrate out fermionic degrees of freedom 



G- 



, where G is the Green-Keldysh 



and arrive at the effective action iS = 2Tt In 
function for our system. Expanding the action in powers of the tunnel Hamiltonian we 
obtain S = Sl,r + Sf^^ + S'P'' + where the term Sl,r. describes the action of the left 
and right conductors, Si^^ is Ambegaokar-Eckern-Schon action [7] and 



iSi 



(2) 



J dti . . . dti J dxi . . . dx4 



i,j,k,l=F,B' 

X GL,j{X,;X2){-iye-^^^^'M^,) X Gn^.kiX^, XsK-lfe'^'^^'Mxs 
X GL,ki{X3;X4){-iye-''^^^'M^4) X Gn,u{X4;X,){-l)Y'^^^'Mxi(}i 
Here X = (t,x), ^Pf{b) is the phase variable on the forward (backward) branch of the 
Keldysh contour, Gr are 2x2 matrix Green-Keldysh functions in the left and right 
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conductors (r = L,R) and we use the convention (—1)'^ = —1, (—1)^ = 1- We assume 
that Gl,r have the equihbrium form Gr = Gf^Fi — F2G^, where G^'^ are retarded and 
advanced Green functions, 



hiE) -- 



f{E) is the Fermi function and h{E) 

Our next step amounts to averaging the products of retarded and advanced 
propagators in the action (j4]) over disorder in each conductor separately. We have 
(see, e.g., [8]) 

(Gf(Xi,X2)Gf(X3,X4)) = (Gf(Xi,X2))(Gf(X3,X4)) 



h{E) - 




h{E) - 


m 




fiE) 


-h{E) - 


~KE) 


1 - f{E) 





+ 27rXiw(|ri - r^\)w{\r2 - r^\) Vl [ti - ta; 
x5(ti - t2 + - ^4) + 2'KNLw{\ri - r^i\)w{\r2 - 



X Cl Ul - t2\ 



6{ti-t2+t3-t4), 



(5) 



2 ' 2 

where N^, V^it, r, r') and Ciit, r, r') are respectively the density of states, the diffuson 
and the Cooperon in the left conductor, w{r) = e~^^'^'''' sin kpr/kpr, kp and Ig are 
respectively the Fermi wave vector and elastic mean free path. The same averaging 
procedure applies to the right conductor. 

Finally we assume that the transmission amplitude t{x) is random, quickly 
oscillating real function. Averaging over these oscillations yields t{x)t{y) = 6{x — 
y)gt{x) I'&'k'^NlNji, where gt{x) is the local conductance of the barrier. After all these 
steps Eq. (jl]) reduces to a sum of different terms. Here we will select only the terms 
responsible for weak localization which involve the product of two Cooperons Cl and 
Cr. Collecting all such contributions we obtain 



iSwL = -i 
9t{x)gt{y) 



X 



j dti . . . dt^ J dTidT2 j dxdy 

Ciiti -Ti,y, x)CR{t2 - Ta, X, y) x e'*^*^ 



M) 



sm 



X 



. i^p (to) . tX) (£9) 

- t2)e-'^ + /(n - 12)^^ 

h{T2 - t3)e'^f{ti + t3 - t4 - Ti) 



-/(r2-t3)e 

_ .■£-{t4) 



'^h{ti + t3-t^-n) 



X 



/(-ti +t2 + ti- T2) + e 



H-ti +t2+t^- T2) 



(6) 



Weak localization in quantum dots 5 

Here we defined "classical" ip~^ = (v?f + V^b)/2 and "quantum" (p~ = (fp — ipB phases 
and introduced 

t^) = - ip+it2) + ip+its) - ^+(t4) 

and f{t) = J{dE/2n) f{E)e-'^^ = 6{t) - h{t). The action (El) fully accounts for the 
effects of electron-electron interactions on WL via the fluctuating phases ip^. 

In order to find the WL correction to the current across the central barrier Iwl we 
make use of the following general formula 



iwL = «e 



J 6lp 



5(p~ 

In the limit ql^r ^ l,gt this integral remains Gaussian in ip^ at all relevant energies and 
can easily be performed. The effective expansion parameter in this case is g}/ gLdR 1. 
Combining Eqs. ([6]) and ^ and introducing the average voltage at the barrier V we 
find 

Iwl = o Re / dxdy gt{x)gt{y) dEduJiduj2duj3 

X CR{-uj2,x,y)CL{-uJ3,y,x) X h{E -uj2)f{E + eV + UJ3-UJ1) 
X [f{E + eV - LUi)h{E)P^{iu,,U2,LU3) + f{E + eV - u^)f{E)P2{u^,uj2.uj3) 
+ h{E + eV - uJi)h{E)P2{uJi, UJ3, UJ2) 

+ h{E + eV- uj,)f{E)P3{iOi,uj2, 003)] - {V^ -V}. (8) 

Here Cl,r and Pj {j = 1, 2, 3) are the Fourier transforms of respectively the Cooperons 
Chfiif) and the functions 

Vj{tiMM) = ^M-^i^iMMWAtiMM), (9) 

where 

T = F{t, + tg) + F{t3) + F{t^ + t2) + F{t2) 

-F{ti+t2+t3)-F{t2-t3) (10) 

and F{t) = {{0{t) — (^(0))^)/2 coincides with the phase correlation function ([2]). The 
terms Qj read 

= e-i[K(t2)+K{t3)+K{\t2-t3\)] 
X {2e'[^(l*l+*2+t3l)+^(tl+t3)+-f^{tl+<2)] _ ^i[K{ti+t2+t3)+K{\ti+t3\)+K{\ti+t2\)]^ 
Q2 = e^i^(\^l+'t2+t3\)-K{t2)-K{\t3\)]^i[K{ti+t3)-K{\ti+t2\)-K{t3-t2)] 
23 = Qill<(^l+-t2+t3)-K{\t2\)~Ki\t3\)]^-i[K{ti+t3)+K{ti+t2)-K{\t3-t2\)] 

where K{t) = i ([<^(0), 0(t)]) is the response function. Eqs. fl8l)- f|TTl) represent the central 
result of our paper. They fully determine WL correction to the current in our system. 
The non-interacting result is reproduced by the first two lines of Eq. (jH]) before the 
square brackets, while the terms in the square brackets exactly account for the effect of 
interactions. The same result follows from the non-linear a- model approach [9]. 
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Our result demonstrates that the whole effect of electron-electron interactions 
is encoded in two different correlators of fluctuating phases F{t) and K{t). These 
correlation functions are well familiar from the so-called P(i?)-theory [3 [ID]- They read 

^^^^ = ^ y ^o^coth — Re[ZM] , (12) 

:r^/ X 9 f duj ^ ,,sinu;t 

K{t) = ^ \ — Re[Z(^)] , (13) 

where Z{lS) is an effective impedance "seen" by the central barrier. Both functions 
f lT^ and f[T^ are purely real and, hence, \Qj\ < 1. At times r^c < 1^1 < iZ-^Th (an 
effective RC-time trc will be defined later) we obtain F{t) ~ ^ (in + 7) and 

K{t) ~ ^signt, where gz = 27i/e'^Z{0) = go + gt, g^^ = gl^ + gR^ + e^Rs/^ix and 
7 ^ 0.577 is the Euler constant. We observe that while F{t) grows with time at any 
temperature including T = 0, the function K{t) always remains small in the limit gz ^ ^ 
considered here. Hence, the combination (fTOjl should be fully kept in the exponent of 
^ while the correlator K{t) can be safely ignored in the leading order in 1/ gz- Then 
all = 1, the Fermi function f{E) drops out from the result and we get Iwl = GwlV 
where 

2 

Sn-^NLNji J J J 

X 9t{x)gt{y)CL{t2,y,x)CR{ts,x,y) 

Identifying ^2 = tf ~ t and t^ = t — ti we observe that the exponent in the third line of 
Eq. (fT^ exactly coincides with the expression ([1]) derived from simple considerations 
involving electrons propagating along time-reversed paths in an external fluctuating 
field. Thus, in the leading order in l/gz the WL correction Gwl is affected by electron- 
electron interactions via dephasing produced only by the "classical" component ip^ of 
the fluctuating field which mediates such interactions. Fluctuations of the "quantum" 
field (p~ turn out to be irrelevant for dephasing and may only cause a (weak) Coulomb 
blockade correction to be considered below. 

It is worthwhile to point out that a similar conclusion was previously reached for 
spatially extended disordered conductors within a different approach [11] . We also note 
that a close relation between the results [H] and the P(i?)-theory [3, [10] was already 
demonstrated earlier [12]. Our present results make this relation even more transparent. 

Our further calculation is concentrated on a system of two (identical) dots depicted 
in Fig. 1. For simplicity the outer barriers are supposed to be open, g^ ^ = g ^ 1, gt 
and Rs — > 0. Then the Cooperons take a simple form CL,R{t) = e~*'''^"/V, where V 
and td are respectively the dot volume and dwell time. We also define the effective 
impedance seen by the central tunnel junction 

/ X 47r 1 / Td 1 1 \ , , 

Z(u)=i- [— r + , 15 
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Figure 2. Temperature dependence of WL correction Gwl (a) and dephasing time 
Tcp (b) for td/trc = 100. 



with the real part 



An 

ReZ{u) = — 



T 



+ 



^2 



(16) 



where 1/r = I/td + Tnc = T^/gEc, Ec = 6^2(0 + Cg + 2Cj) and C, Cj and 

Cg are the capacitances of respectively left (right) barriers, the central junction and the 
gate electrode. Substituting the Cooperons CL^Rit) and the correlator F{t) (fT2|) . (fT6ll 
into Eq. (fn|) we observe that contribution of 5{ijj) in Eq. (ITB]) drops out. Performing 
the time integrals we arrive at the final expression for the WL correction Gwl {T) in the 
presence of electron-electron interactions: 



G 



WL 



dt2dt3 e~(*2+t3)/TI3g~2F(i2)-2F{t3)+F(i2+t3)+F(t2-i3)^ (j^^) 



where S is the dot mean level spacing. This result is plotted in Fig. 2a 
demonstrating that interactions suppress Gwl(T) below its non- interacting value [S] 
G^Sl = -2e'9?/7rg\ 

Let us define u = td/trc = AEc/S and consider the limit of metalhc dots m ^ 1. 
At Tto < 1 the WL correction saturates to 

--(0) I ~ -(2/uf', g> 



G 



WL/ \Gwl\ ^ 

\G'Sl\ - 



Gwl/\G'^l\- -^?/2m, 1<^?<8, 
whereas aX g / tjj < T < \ / trq and for (7 > 8 we find 



G 



WL 



\G 



(0) I 

wl\ 



9__^ \j2nM^ 
4 7 (7rTri3)i-8/5- 



(19) 



Let us phenomenologically define the electron decoherence time by taking the 
Cooperons in the form Ci^Rit) = e~*/^^-*/""^/V which yields H] Gwl/G^^^ = (1 + 
td/t^)^'^. Resolving this equation for we obtain 



r,. 



WL 



(20) 



which yields = gTD/A\Yi{2Ec/5) for Tr^ < 1 and ^ > 81n(M/2). Eqs. ([I8D-([20D, 
although not directly applicable to a single quantum dot, account for key features of 
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the dependence r<^(T) (Fig. 2b) observed in various UCF experiments [ISIEIIIS] with 
quantum dots |T6]. At higher temperatures we find r^p oc with non- universal g- 
dependent power z/ < 1/2, while at lower T < 1/td the electron decoherence time 
saturates to a constant in agreement with the observations [ISj [m [15] . It was pointed 
out [15] that the available experimental values of r^(0) scale as r<^(0) ~ td for a variety 
of dot sizes and dwell times td varying by ~ 3 decades. Our result fll8|20p should be 
consistent with this scaling provided at low T the right-hand side of Eq. (I20l) remains 
of order one. 

Note that the phenomenological definition of is identical to that used before 
in Ref. [9] where we also demonstrated that for an arbitrary array of quantum 
dots our expression for the weak localization correction determines the system 
magnetoconductance if we substitute 1/r^ — >• l/r^, + 1/th, where th oc I/H"^ is the 
electron dephasing time due to the external magnetic field H. Thus, our definition of r<^ 
is fully consistent with the standard procedure of extracting the electron dephasing time 
from the magnetoconductance curves. Furthermore, it is straightforward to demonstrate 
[9] that, e.g., in the case of quasi-ld arrays of quantum dots our definition for just 
yields the standard result for the magnetoconductance of a diffusive wire, cf. Eq. (60) 
of Ref. [9]. 

We also would like to emphasize that there exists no contradiction between the 
definition of r^p adopted here and the fact that no dephasing occurs for electon paths 
confined within a single quantum dot, as discussed in the beginning of our manuscript. 
As it was demonstrated, electron dephasing occurs as soon as time-reversed paths cross 
the central barrier twice and return to the initial point inside the dot (see Fig. 1). In the 
presence of fiuctuating electromagnetic potentials (dropping across the central barrier) 
the forward path and its time-reversed counterpart pick up different random phases. 
After averaging over both fluctuating flelds and electron paths one arrives at a decaying 
in time contribution to the Cooperons which is just captured by our phenomenological 
definition. Of course, other definitions of can also be employed. However, our basic 
conclusion about non-vanishing electron dephasing by electron-electron interactions 
down to T — > will not be sensitive to any particular definition of r<^, since this 
conclusion is based on the result ( fT4l) demonstrating the interaction-induced suppression 
of the WL correction to conductance (as well as of the magnetoconductance, cf. Ref. 
[9]) at any temperature including T = 0. The basic physics behind this result is exactly 
the same as that already elucidated by the well known P(i?)-theory [TJ ^0\: tunneling 
electrons can exchange energies with an effective electromagnetic environment. This 
process results in broadening of the distribution function for such electrons even at 
T = which inevitably yields electron dephasing. 

Finally, it is instructive to establish the relation to the ordinary perturbation theory 
in the interaction which is reproduced by formally expanding our exact result ([8]) to the 
first order in Z{u)). We obtain 

IwL = GI^\V + 6I^^{V) + 6I^dV), (21) 
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where 



crF (f^lS^ f duj ^R^Z{uj) uj 

^^^^ = -^'^J 2^^ 2T 



(22) 




+ ilmZ{uj)CLiuj)CRiuj)}. 



(23) 



Here we defined the function 



W = {uj + eV) coth 



uj + eV 



{u — eV) coth 



u-eV 
2T 



2T 



and Fourier transformed Cooperons Cl{uj) = Cr{u!) = to/{1 — iujtd)- The two terms 
SI^j^ and are hnear in respectively F{t) and K{t). 

Exactly the same results ( l2T|) -(l23l) are reproduced from the first order diagrammatic 
perturbation theory in the interaction. In order to observe the equivalence of the two 
approaches one should keep in mind that F{t) is proportional to the Keldysh component 
of the photon Green function, while K{t) is proportional to the retarded photon Green 
function. One should also remember that the photon Green function in our model 
is coordinate independent in both quantum dots. One can actually demonstrate that 
the terms oc Gl{u!)Gr{0), Gl{0)Gr{u!) come from the so-called "self-energy" diagrams, 
while the terms oc Gl{—uj)Gr{uj) emerge from the "vertex" diagrams. 

The term 6I^r^ represents the Coulomb blockade correction to and is entirely 
different from the dephasing term In contrast to the latter, the term SI^j^ is 

non-linear in V describing the standard Coulomb offset at large V and turning into 



for Tr^ > 1 in the linear in V regime. Thus, the Coulomb blockade correction remains 
small [IT] in the metallic limit g ^ 1. We also note that SI^rl involves the combination 
l — 2f{E) = tanh(£'/2T) which enters only in the first order in the interaction. As in the 
case of spatially extended conductors, at T = some terms contained in 5/^^ partially 
cancel similar contributions to 51^ i^. This cancellation, however, remains incomplete 
and, as demonstrated by our exact result, by no means implies the absence of electron 
dephasing at T — > 0. More information on the debates on low temperature decoherence 
by electron-electron interactions can be obtained, e.g., from Refs. [18] and further 
references therein. Without going into details, we would only like to emphasize that our 
present manuscript does not make any use of the techniques introduced in our previous 
works on decoherence in disordered conductors and, hence, is formally independent on 
those. 

In summary, we have non-perturbatively treated the effect of electron-electron 
interactions on weak localization in relatively short metallic conductors. The most 




(24) 
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significant effect of interactions is electron decoherence which persists down to T = 
and - in agreement with experiments [131 El [IS] - yields saturation of r<^ at T < 1/r/). 
The physics behind this effect is exactly the same as that discussed, e.g., within the 
well known P(i?)-theory [3 [10]. It is also worth pointing out that very recently [19] we 
generalized our present approach to arbitrary arrays of quantum dots and derived the 
expression for t^q which describes both weakly and strongly disordered conductors and 
quantitatively explains numerous experimental data available to date. In the case of 
weakly disordered conductors our results [I9] match with those derived previously [H] 
by means of a different technique. 

This work was supported in part by the EU Framework Programme NMP4- 
CT-2003-505457 ULTRA-ID "Experimental and theoretical investigation of electron 
transport in ultra-narrow 1-dimensional nanostructures" and by RFBR Grant 06-02- 
17459. 
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